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Super-Kamiokande is a 50 kiloton water Cherenkov detector located at the Kamioka Observatory 
of the Institute for Cosmic Ray Research, University of Tokyo. It was designed to study neutrino 
oscillations and carry out searches for the decay of the nucleon. The Super-Kamiokande experiment 
began in 1996 and in the ensuing decade of running has produced extremely important results in the 
fields of atmospheric and solar neutrino oscillations, along with setting stringent limits on the decay 
of the nucleon and the existence of dark matter and astrophysical sources of neutrinos. Perhaps 
most crucially, Super-Kamiokande for the first time definitively showed that neutrinos have mass 
and undergo flavor oscillations. This chapter will summarize the published scientific output of the 
experiment with a particular emphasis on the atmospheric neutrino results. 



I. INTRODUCTION AND PHYSICS GOALS 

The Super-K collaboration is the combination of two 
previous successful collaborations. The first was the 
Kamiokande [j], 0] experiment which was located in the 
same mine as Super-K and had a fiducial mass approxi- 
mately 20 times smaller than Super-K. The second was 
the 1MB experiment [H, 0] which was located in the Mor- 
ton Salt mine in Ohio. Grand Unified models such as 
SU(5) [5j predicted that the proton would decay at a 
rate visible by modest size detectors and both of these 
experiments were originally built to search nucleon decay 
into the mode favored by SU(5) which is p — > e + 7r°. 

Although neither of these experiments observed the de- 
cay of the proton, they did measure a statistically signif- 
icant anomaly in the expected background to the proton 
decay search from neutrino interactions on the water in 
the tanks. One explanation for this effect was that some 
of the neutrinos were oscillating into an un-observable fla- 
vor and thus giving less background than expected. At 
the same time two detectors made of iron the NUSEX [6] 
and Frejus 01 experiments saw a result which was con- 
sistent with the expectation but with much lower statis- 
tics. Super-Kamiokande was designed to definitively de- 
termine whether or not oscillations were indeed taking 
place. 

Additionally, by scaling up the size from previous de- 
tectors, Super-Kamiokande offered new hope to finally 
observe the decay of the nucleon and also to try to answer 
the crucial question of whether neutrinos produced in the 
nuclear burning in the sun oscillated into non-detectable 
flavors on their way to the earth. Previous generations 
of experiments had not seen as many neutrinos from the 
sun as predicted by solar models. With a large mass, 
good energy resolution, and the ability to point neutrinos 
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back to the sun in real-time, Super-Kamiokande was de- 
signed first of all to confirm the effect with high statistics 
and then to determine what the parameters of oscillation 
were. 



II. THE SUPER-KAMIOKANDE DETECTOR 

Super-Kamiokande is a 50 kiloton water Cherenkov de- 
tector located at the Kamioka Observatory of the Insti- 
tute for Cosmic Ray Research, University of Tokyo. Fig- 
ure[TJshows a diagram of the Super-Kamiokande detector. 
The detector is in the Mozumi mine of the Kamioka Min- 
ing Company in Gifu prefecture, in the Japanese alps. 
Super-K consists of two concentric, optically separated 
water Cherenkov detectors contained in a stainless steel 
tank 42 meters high and 39.3 meters in diameter, holding 
a total mass of 50,000 tons of water. The inner detector 
is comprised of 11,146 Hamamatsu R3600 50 cm diam- 
eter photomultiplier tubes, viewing a cylindrical volume 
of pure water 16.9 m in radius and 36.2 m high. 




FIG. 1: An overview of the Super-Kamiokande detector site, 
under Mt. Ikeno from Ref. @]. The cutaway shows the inside 
lined with photomultiplier tubes comprising a photo-cathode 
coverage of about 40%. 
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As described more fully in Ref. ||, the detector is cal- 
ibrated in energy at the 2% level over enegies ranging 
from the MeV to the tens of GeV range. This careful 
calibration is key to the successful extraction of physics. 
It relies both on natural calibration sources such as the 
expected energy deposit of muons created by the cos- 
mic rays and the decay of neutral pizeros produced by 
neutrino interactions inside the Super-Kamiokande tank, 
and on artificial sources including lasers, a Xenon light 
source, a linear accelerator @ and a 16 TV source [icj . 
Low energy radioactive background were also carefully 
studied, a description of the measurement of radon con- 
centrations at Super-Kamiokande is given in Ref. [111 ]. 

Data from the detector is first collected by an on-line 
data acquisition system and then, after a calibration step, 
passed into several streams of reduction for the various 
analyses. Selection steps are performed to remove the 
background from non-neutrino induced interactions. For 
example, the outer detector region is used as a veto to 
reject cosmic-ray muons. Using the time and charge in- 
formation at each photo-tube, reconstruction algorithms 
are applied to the data to determine a vertex for the 
Cherenkov light and any rings associated with the parti- 
cles of the interaction. Additional likelihood-based algo- 
rithms are used to determine the properties of the par- 
ticles that generated the light including their type, mo- 
menta and directions. These reconstructed physics quan- 
tities are then used for analysis. More detail about the 
acquisition and analysis reconstruction techniques can be 
found in Refs. HOj and [H|. 

The Super-Kamiokande running periods are divided 
into three parts. The first, SK-I ran from 1996 to 2001. 
In November of 2001 an accident destroyed much of the 
photo-tubes of Super-K and, after a year of rebuilding the 
detector with half of the previous photo-tube density, SK- 
II ran for approximately 800 days. In June of 2006 after 
restoring Super-Kamiokande to its full photo-tube den- 
sity a period of running known as SK-III began. During 
the SK-I and SK-II running period Super-Kamiokande 
acted as the target for the long-baseline K2K experi- 
ment. Starting in 2009, Super-Kamiokande will once 
again be the target of an accelerator produced neutrino 
beam when the T2K experiment begins. 



In the sections that follow, papers from all of these sub- 
jects will be reviewed. Particular attention will be payed 
to the history of the published papers in the atmospheric 
neutrino analysis. 



A. Atmospheric Neutrino Oscillations 
The sub-GeV R ratio Q 

In the previous results from the 1MB and Kamiokande 
experiment it was observed that the flavor ratio of neu- 
trinos below 1 GeV did not agree with expectation. If 
one took the results at face value, either there were more 
electron neutrinos then expected, or too few muon neu- 
trinos. In order to study the question experimentally a 
measurement was made of: 



R = (p / e) D AT A / (a*/ s)mCi 



(1) 



which served to cancel uncertainties in neutrino flux and 
cross-sections. In Eqn. [TJ (n/e) means the ratio of the 
number of measured neutrino interactions which are in- 
ferred to have come from electron and muon neutrinos 
respectively. The ratio is calculated separately for the re- 
constructed data and MC. If there is perfect agreement 
between data and expectation the expected value of R 
is one. In Super-Kamiokande, the measured value of R 
was: 



R = 0.61 ± 0m(stat.) ± 0.05(sys.). 



(2) 



This was a statistically significant result and the col- 
laboration concluded in Ref. 



14|: 



"The first measurements of atmospheric neutrinos in 
the Super-Kamiokande experiment have confirmed the 
existence of a smaller atmospheric v^jv e ratio than pre- 
dicted. We obtained R = 0.61 ± 0.03(stat.) ± 0.05(sys.) 
for events in the sub- GeV range. The Super-Kamiokande 
detector has much greater fiducial mass and sensitivity 
than prior experiments. Given the relative certainty in 
this result, statistical fluctuations can no longer explain 
the deviation of R from unity. " 



III. PUBLISHED RESULTS FROM 
SUPER-KAMIOKANDE 

The published scientific output of the Super- 
Kamiokande experiment can be roughly divided into four 
main categories: 

1. Studies of atmospheric neutrino oscillations 

2. Studies of solar neutrino oscillations 

3. Searches for the decay of the nucleon 

4. Searches for astrophysical sources of neutrinos 



The multi-GeV R ratio [IE 

The previous result relied on events which had visi- 
ble energy less that 1.33 GeV deposited in the Super-K 
tank (so called sub-GeV events). Although less numer- 
ous, multi-GeV events were also expected to oscillate and 
had the extra advantage that at high-energy the outgo- 
ing lepton direction closely followed the incoming neu- 
trino direction. Since the neutrino oscillation probabil- 
ity is a function of both the distance traveled and the 
energy, knowing the direction of the incoming neutrino 
allowed the separation of neutrinos into bins of angular 
zenith. Based on the results of the previous experiments 
it was expected that neutrinos of a few GeV would need 
to travel thousands of kilometers before they oscillated. 
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Also, if muon neutrino were oscillating into tau neutri- 
nos, the vast majority of them would not have the energy 
necessary to interact and produce a tau lepton. In this 
case, the expectation was that there would be a deficit of 
muon interactions coming from below. 

From the data analyzed, the R ratio for the multi-GeV 
events was reported in Ref. [l5[ to be: 



R = 0.66 ± 0M(stat.) ± 0.08(sys. 



(3) 



confirming the result seen in the sub-GeV sample. 
Crucially, the expected zenith suppression was also seen. 



Evidence for oscillations [If 

In 1998, the paper "Evidence for oscillation of atmo- 
spheric neutrinos" was published 16]. In this paper, 535 
days of data were analyzed and subdivided into sub-GeV 
and multi-GeV e-like and mu-like events. In this data set, 
a strong suppression of upward going event from muon- 
neutrino interactions was observed. This was quantified 
in terms of an asymmetry defined as A = (U—D)/ (U+D) 
where U is the number of upward-going events (—1 < 
cos0 < —0.2) and D is the number of downward-going 
events (0.2 < cos0 < 1). Due to the isotropic nature of 
the cosmic rays this asymmetry is expected to be close 
to zero. For the electron-neutrino sample this was found 
to be the case. For the muon-like events the asymmetry 
was found to be: 



A = 



(U-D) 
(U + D) 



-0.296 ± 0.048(siai.) ± 0.01(sys.) (4) 



which deviates from zero by more than 6 standard devi- 
ations. Figure [2] taken from Ref. [l6| shows this asym- 
metry as a function of momentum for both the e-like and 
mu-like samples. The deficit of high energy upward-going 
mu-like events is clearly seen. 

A fit was performed over the oscillation parameter 
space to the Vn «-> v r oscillation hypothesis. There were 
eight systematic uncertainty terms that were allowed to 
vary within their known ranges. Fits were also performed 
to the v n <-* v e hypothesis but the data did not fit this 
hypothesis well. The data-samples and the results of the 
fit are shown in the Fig. [3] also from Ref. [16fl . 

This analysis resulted in the first allowed region for at- 
mospheric neutrino oscillations from Super-Kamiokande, 
resulting in a somewhat lower allowed region then was 
previously obtained from the Kamiokande collaboration. 
In this paper the conclusion was: 

Both the zenith angle distribution of /i-like events and 
the value of R observed in this experiment significantly 
differ from the best predictions in the absence of neutrino 
oscillations. While uncertainties in the flux prediction, 
cross sections, and experimental biases are ruled out as 
explanations of the observations, the present data are 
in good agreement with two-flavor <-> v T oscillations 
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FIG. 2: The (U - D)/(U + D) asymmetry as a function of 
momentum for FC e-like and /i-like events and PC events. 
While it is not possible to assign a momentum to a PC event, 
the PC sample is estimated to have a mean neutrino energy 
of 15 GeV. The Monte Carlo expectation without neutrino 
oscillations is shown in the hatched region with statistical 
and systematic errors added in quadrature. The dashed line 
for /i-like is the expectation for «-» v T oscillations with 
(sin 2 26 =1.0, Am 2 = 2.2 xlCT 3 eV 2 ). Figure and caption 
from Ref. Q3]. 



at 90% confidence level. We conclude that the present 
data give evidence for neutrino oscillations. 

It should be noted that the allowed region found in 
1998 is consistent with all of the later more refined 
analysis, which have also all been consistent with each 
other. 



with sin" 2 26 > 0.82 and 5 x 10" 4 <Am i < 6 x lfT 



iV 2 



Upward-going [17] and stopping [18] muon sam- 
ples 

Also analyzed separately in Refs [13] and [l8| were up- 
ward through-going and stopping muon samples. Muons 
are created by neutrino interactions under the Super-K 
tank and then travel upwards, either passing through or 
stopping inside the detector. The through-going events 
were of a higher energy than the events analyzed in 
Ref. [l6| which were completely contained inside the 
Super-K tank, and the stopping events were of compara- 
ble energy to the events previously considered in which 
the produced muon escaped the tank. The reduction pro- 
cess for these events was somewhat more involved since it 
was more difficult to remove backgrounds from entering 
downward-going muons. 

In these two papers, the two sub-samples were fit to 
the oscillation hypothesis as was their ratio as a function 
of zenith angle. It was found that they also gave results 
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FIG. 3: Zenith angle distributions of /^-like and e-like events for sub-GeV and multi-GeV data sets. Upward-going particles 
have cosO < and downward-going particles have cosO > 0. Sub-GeV data are shown separately for p < 400 MeV/c and 
p > 400 MeV/c. Multi-GeV e-like distributions are shown for p < 2.5 GeV/c and p > 2.5 GeV/c and the multi-GeV p,-like are 
shown separately for FC and PC events. The hatched region shows the Monte Carlo expectation for no oscillations normalized 
to the data live-time with statistical errors. The bold line is the best-fit expectation for <-> v T oscillations with the overall 
flux normalization fitted as a free parameter. Figure and caption from Ref. [in- 



consistent with neutrino oscillations as presented in 
Ref. 0. 

The East- West [l9| effect 

In Ref. [l!| the east-west effect was seen in neutrinos 
for the first time. This effect, first since in cosmic ray 
showers, is a direct confirmation of the mostly positively 
charged nature of the cosmic rays. Low-energy cosmic 
rays are deflected by the Earth's magnetic field distorting 
the measured azimuthal spectra. By selecting a subset of 
events sensitive to the Earth's magnetic field this effect 
was observed as shown in Fig. 2] taken from [l^| . 

While the previous analyses explored only zenith 
angle distortions caused by oscillations, by confirming 
the expected azimuthal distortion it was demonstrated 
that the geomagnetic effects on the production of GeV 
energy neutrinos was well understood. 

Tests for non-standard oscillation models [20] 

As more data was collected and the techniques were 
improved much effort was put into trying to test if any 
other explanation was as good as the «-» v r at explain- 
ing the observed data. 

In Ref. [13] with over 1000 days of data the fully and 
partially contained events were fit together along with 
the upward-going muons and a multi-ring sample which 
had been enriched with neutral-current interactions. 
In «-> v r oscillations, a simple disappearance of the 
charged current events would be expected since tau 



neutrinos still have neutral current interactions. On the 
other-hand, if the oscillations were to sterile neutrinos 
which had no interactions at all, two additional effects 
would arise. First, there would also be a reduction 
in the neutral current sample. Second, because of a 
difference in the forward scattering cross-sections in 
the two cases, the normal charge current oscillations 
would be suppressed at high energies. This suppression 
is known as the "matter effect." A fit was done to 
the data samples looking for these effects and was 
not seen. Instead, the data fit the standard <-» v T 
hypothesis extremely well and excluded oscillations into 
sterile neutrinos at greater than the 99% hypothesis level. 

Observation of the oscillation pattern [21] 

In 2004 in Ref. [2l[ the standard oscillation hypothesis 
was given strong confirmation by the observation of a 
sinusoidal pattern in the oscillation pattern. 

Since the oscillation equation contains a sin function, 
in principle one should be able to see a sinusoidal dip 
if the data is plotted as a function of L/E. However, in 
the standard analyses, low resolution in reconstructed 
L (distance from production) and E (neutrino energy) 
washes out the effect if the data is plotted in these 
combinations of variables. In Ref. [2l[ a subset of events 
were selected which had high resolution in these vari- 
ables and after plotting the observed suppression relative 
to the expectation as a function of L/E the tell-tale 
dip from oscillations was observed. Additionally, by 
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FIG. 4: Azimuthal angle distributions of e-like and /x-like 
events. The crosses represent the data points, the histogram 
drawn by solid line (dashed-line) shows the prediction of the 
Monte Carlo based on the flux of [two different flux mod- 
els]). Data are shown with statistical errors. The Monte Carlo 
has f times more statistic than data. The Monte Carlo his- 
togram is normalized to the total number of the real data. <j> 
represents the azimuthal angle. (6 = 0, 7r/2, n and 3n/2 shows 
particles going to north, west, south, and east, respectively. 
Figure and caption from Ref. [r9T ]. 
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FIG. 5: Ratio of the data to the MC events without neutrino 
oscillation (points) as a function of the reconstructed L/E 
together with the best-fit expectation for 2-flavor «-» v T 
oscillations (solid line). The error bars are statistical only. 
Also shown are the best-fit expectation for neutrino decay 
(dashed line) and neutrino decoherence (dotted line). Figure 
and caption from Ref. [2lj ] 



significance. Figure \5\ demonstrating the effect is taken 
from Ref. 

Measurement of atmospheric neutrino parame- 
ters (3 

In 2005 all of the data from SK-I was refitted. This 
time the fit was done jointly with all of the data includ- 
ing upward-going and multiple ring-samples together. 
Ref. [12j describes this analysis in detail along with de- 
tails of the atmospheric neutrino Monte Carlo and recon- 
struction algorithms and reduction, 
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FIG. 6: The 68, 90 and 99 % confidence level allowed oscilla- 
tion parameter regions obtained by an L/E analysis [2lJ and 
by the 2005 complete data-set analysis are compared. Note 
the use of the linear y-scale in the figure. Figure from Ref [l2| ]. 



Thirty-nine systematic uncertainties were accounted 



for in the fit, and the best fit for v u 



oscillations 



with sin 2 26> = 1.00 and Am 2 = 2.1 x 10~ 3 eV 2 was 
found. Figure [Jj taken from Ref. [12| and shown on 
the next page, displays all of the data samples and the 
results of the fit. Over 15,000 atmospheric neutrino 
events were used in the analysis. In Fig. [S] the final 
allowed region is shown along with the allowed region 
that is found by the L/E analysis from Ref. [20 ] . Both 
of these analyses give consistent results. The better 
constraint in sin 2 29 in the zenith angle analysis is 
due to higher statistics in the sample, while the better 
constraint in Am 2 in the L/E analysis is due to finer 
binning in the fit. Later, as yet unpublished analyses 
combine the best feature of both techniques to extract 
the maximum information. 



comparing the x 2 of the data to the standard oscillation 
scenario to that of neutrino decay and decoherence those 
other exotic hypothesis could be rejected with high 



Three flavor oscillations [22] 



In Ref. [12| the data was fit assuming that there were 
only oscillations between two flavors of neutrinos. How- 
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FIG. 7: The zenith angle distribution for fully-contained 1-ring events, multi-ring events, partially-contained events and upward 
muons. The points show the data, box histograms show the non-oscillated Monte Carlo events and the lines show the best-fit 
expectations for <-» v T oscillations with sin 26» = 1.00 and Am 2 = 2.1 x 10 3 eV 2 . The best-fit expectation is corrected by 
the 39 systematic error terms, while the correction is not made for the non-oscillated Monte Carlo events. The height of the 
boxes shows the statistical error of the Monte Carlo. Figure and caption from Ref. [12J. 



ever, we know there are three flavors of neutrinos and 
one can write down a three-flavor mixing matrix with 
three mixing angles. The two that have been measured 
(#12 and #23) control solar and atmospheric mixing re- 
spectively. The third, as yet unmeasured, mixing angle 
813 would allow all three of these neutrinos to oscillate 
into each-other, producing small amounts of high energy 
electron neutrino appearance in the upward going data. 
The question of the value of #13 is particularly pressing 
since if it is extremely small or zero it will not be possible 



to measure CP violation in the neutrino sector using the 
oscillation technique. 

In order to measure or constrain the value of #13 
using the Super-K data set, a fit was performed using 
full three flavor oscillation probabilities, taking into 
account the resonance effects that can occur in the 
earth for such oscillations. The data was binned so as 
to be maximally sensitive to upward going electrons in 
the relevant energy regions. As detailed in Ref. [22j no 
evidence for non-zero #13 was found. A plot showing the 
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Super-K 6*13 allowed region overlaid with the exclusion 
region for the CHOOZ reactor experiment is shown in 
Fig. EI 
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FIG. 8: Allowed regions of #13 as found in the SK-I data set in 
the three flavor oscillation analysis. 90 % (thick line) and 99 % 
(thin line) confidence level allowed regions are shown in Am 2 
vs sin 2 #13. Normal mass hierarchy (Am 2 >0) is assumed. The 
shaded area in the bottom figure shows the region excluded 
by the CHOOZ reactor neutrino experiment. From Ref [22I ]. 



Search for tau lepton appearance [2 

None of the previous analyses explicitly searched for 
the appearance of the tau lepton from oscillations. The 
number of tau leptons in the data sample is expected 
to be small. Only on the order of 80 events with tau 
leptons were expected to be produced in fiducial volume 
of Supcr-Kamiokandc during the SK-I running period. In 
Ref. [231 ] an explicit search for tau lepton appearance was 
undertaken. 

This search relied on the fact that events containing 
heavy tau leptons which decay hadronically, decay 
symmetrically with more pions in the final state then 
the charged current deep-inelastic events which form the 
background to the search. A statistical separation of 
signal to background was made using both a likelihood 
and a neural-network. The background was normalized 
with downward going events by taking advantage of the 
fact that the entire tau signal comes from below. Fig [9] 
from Ref. t 23j shows the results of the fit. This analysis 
measured a tau appearance signal which was consistent 



with that expected from v„ 



oscillations. 



B. Solar Neutrino Oscillations 

The neutrinos produced in the nuclear burning of the 
sun are of lower energy than atmospheric neutrinos. 



Likelihood 




-0.5 





cose 



0.5 



FIG. 9: The zenith angle distributions for the 
likelihood (top) and neural network (bottom) 
analyses. Zenith angle cos(9 = — 1 (cos# = +1) 
indicates upward-going (downward-going). The 
data sample is fitted after tau neutrino event se- 
lection criteria are applied. The solid histogram 
shows the best fit including v T , and the dashed 
histogram shows the backgrounds from atmo- 
spheric neutrinos (v e and z/ M ). An excess of tau- 
like events is observed in upward-going direction 
(shaded area). From Ref. [231 ] 



Super-Kamiokande is sensitive mostly to neutrinos from 
the S B branch of the pp nuclear fusion chain in solar 
burning. At very low energies the experiment is domi- 
nated by radioactive backgrounds. However, above ap- 
proximately 4 MeV the detector can pick-out the scat- 
tering of solar neutrinos off atomic electrons which make 
Cherenkov light in the tank. The S B and rarer HEP 
neutrinos have a spectrum which ends near 20 MeV. 

Unlike previous radio-chemical experiments which re- 
lied on extraction of isotopes, Super-Kamiokande collects 
its data in real time and the electrons which are scattered 
by neutrinos point in a direction that is correlated with 
the sun. Therefore, by plotting the direction between low 
energy events in the tank and the Sun, one can pick out 
a peak containing solar neutrinos. This is clearly demon- 
strated in Fig. [TU1 which the cos of the angle between the 
events and the sun are plotted. 

The first solar neutrino results from Super- 
Kamiokande were presented in 1998 in Ref. [24| 
and reported the measured flux from the first 300 days 
of data. In this paper, only the measured rate was 
reported and Super-Kamiokande confirmed the existence 
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FIG. 10: Angular distribution ol solar neutrino event candi- 
dates. The shaded area indicates the elastic scattering peak. 
The dotted area is the contribution from background events. 
From Ref. [H 



of deficit in neutrinos expected from by the sun. It was 
consistent with the previous Kamiokande measurement 
and reported a flux which was 36% of the standard 
solar model of 1995. For this first analysis, the energy 
threshold was set at 6.5 MeV which was later lowered. 

The next two papers used 504 days of data and 
searched for distinctive signs of neutrino oscillation. In 
Ref (25|, a measurement was made of the solar flux as 
measured at Super- Kamiokande in the day and the night 
separately. During the day, neutrinos from the sun travel 
down from above, traveling through very little earth. 
While during the night they travel from the other side 
of the Earth, in some parts of the year actually pass- 
ing through the outer core. This is relevant because, for 
some regions of oscillation space, a regeneration of the 
original flux can happen inside high density materials. 
In this paper, no statistically significant sign of this was 
seen, thereby excluding parts of the oscillation space pre- 
viously allowed by other experiments. At the same time, 
an analysis of the measured energy spectrum shape was 
published in Ref. [26| using the 504 day data set. It was 
found that a \ 2 analysis gave a a probability of 4.6% of 
being in agreement with the standard solar model. 

In 2001, with more than twice as much data, the anal- 
ysis threshold was lowered to 5 MeV and the systematics 
of the measurement were reduced due to extensive cal- 
ibrations and improvements in the analysis techniques. 
In Ref. [271 ] a precise measurement of the solar flux was 
found. Additionally, the spectrum was found to have no 
significant energy spectrum distortions, and a no statis- 
tically significant day-night effect was seen. The small 
but expected seasonal variation of the flux due to the 
eccentricity of the Earth's orbit was observed and new 
stringent limits on the flux of HEP neutrinos were pre- 
sented. The lack of spectral and zenith angle distortions 
placed strong constraints on the oscillation solutions to 
the solar neutrino problem and in Ref. 28] an oscilla- 



tion analysis was performed using both the spectral and 
flux information from Super-Kamiokande. The result of 
this analysis was a preference for the Large Mixing Angle 
solution. 

In 2002 an analysis including all 1496 days of SK-I was 
published [2!| . The lack of spectral distortion and daily 
variation in the flux strongly constrained the regions al- 
lowed for neutrino oscillation from other experiments, 
leaving only the high-mass LMA and quasi-vacuum re- 
gion allowed. If the Super-K interaction rate and either 
the standard solar model prediction of the a B flux, or the 
rates as measured by SNO were included in the analysis, 
the LMA angle was uniquely allowed at high confidence 
level. This single allowed region could uniquely explain 
the solar neutrino problem. 

The full SK-I data set was also used to search for other 
exotic signals from the Sun. First was the search for anti- 
electron neutrinos in Ref. (30j . The motivation for this 
search was to exclude the possibility that the electron- 
neutrinos from the sun are disappearing do to a spin- 
flavor-precession where neutrinos with a large magnetic 
moment are transformed into their anti-particles in the 
strong magnetic field of the sun. Unlike the standard 
solar neutrino analysis, the reaction here is inverse beta- 
decay off of protons in the nucleus. No signal was found 
and a limit was set of 0.8% of the standard solar neutrino 
flux between 8-20 MeV. 

Next, in Ref. [3l|, searches for periodicity in the solar 
neutrino data using the Lomb test were performed. No 
significant statistical fluctuation was found. Finally, in 
2004, a search for the magnetic moment of the neutrino 
was presented in Ref. [12] • The SK-I data set was used to 
look for energy distortion of electron recoil spectrum. Us- 
ing Super-Kamiokande's allowed oscillation parameters a 
limit of < 3.6 x 10~ 1( Vs was set. If the constraints on 
oscillation parameters from other experiments were also 
included then the limit was reduced to < 1.1 x 10~ 10 /is. 

In 2004, in Ref. 33], the full SK-I data set was once 
again analyzed but this time with a maximum likelihood 
analysis applied to the zenith angle dependence of the 
data. This lowered the statistical uncertainty on the 
day/night ratio by 25% compared to the previous mea- 
surement. This was equivalent to adding three years 
of live-time running if still using the previous method. 
Fig. [11] from Ref. [33| shows the data and best fit spec- 
trum in the LMA region along with the D /N asymmetry 
as a function of energy. 

In 2005 a full detailed description of the SK-I solar 
neutrino analysis was published in Ref. [Hil ]. This paper 
included descriptions of the simulation, reconstruction 
and analysis techniques. In the large mixing angle region 
the day-night asymmetry was found to be: 

_ ($p - <P N ) _ 

- 1.7% ± 1.6%(atat.)tij$(»y*.) ± .04%(Am 2 ) (5) 
which is statistically compatible with zero and is to be 
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FIG. 11: LMA Spectrum (top) and D/N Asymmetry (bot- 
tom). The predictions (solid lines) are for tan 2 f3 = 0.55 



and Am 



6.3 x 10 _b eV 2 with fa, 



0.96 x Standard So- 



lar Model [BP2000] and fa CT = 3.6 x Standard Solar Model. 
Each energy bin is fit independently to the rate (top) and 
the day/night asymmetry (bottom). The gray bands are the 
±1<7 ranges corresponding to the fitted value over the entire 
range 5-20 MeV: A = —1.8 ± 1.6%. Figure and caption from 
Ref. 11. 



compared with the expected asymmetry which ranges 
from -1.7% to 1.0%. 



C. The Search for Proton Decay 

Super-K has also set important limits on the decay of 
the nucleon. The limits from Super-K have ruled out 
SU(5) and the minimal super-symmetric model. The 
published limits on p ~ + e + ir° as predicted in SU(5) are 
found in Ref. [p. 

Many models of super-symmetry predict that the nu- 
cleon should decay with strange quarks in the final state. 
The first 535 days of SK data were analyzed to search 
for p — > DK + and the limits were presented in Ref. [35| . 
In 2005 the full SK-I data set was employed to search 
for SUSY mediated decays of the proton. In the anal- 
ysis presented in Ref. [361 ] stringent limits were set on 



p — > vK + ,n— > v K" , p — > /i ' A " and p - 

A related search for SUSY generated physics is the 
search for neutral Q-Balls as presented in Ref. [37|- Q- 
balls are topological solitons predicted in some SUSY 
models and interact as they pass through the Super-K 
tank leaving a trail of pions in their wake. A search was 
carried out using 542 days of the SK-II data set. No 
evidence for Q-balls were found, and Super-K finds the 
lowest limits in the world for Q-ball cross-sections below 
200 mb. 



i + K° and 



K modes. 



D. The Search for Astrophysical Phenomenon 

Super-Kamiokande has been used to search for sev- 
eral sources of neutrinos from outside of our solar sys- 
tem. Gamma Ray Bursters are among the most lumi- 
nous sources that have ever been observed in the uni- 
verse. Super-K performed a search for neutrinos that 
were in coincidence with the BATSE detector located on 
NASA's Compton Gamma Ray Observatory. The en- 
tire Super-K data sample from 7 MeV to ~ 100 TeV was 
compared with the BATSE online catalog in the period of 
April 1996 to May of 2000. The results were presented in 
Ref. [38| and no statistically significant signal was found. 

A supernova which occurred in the center of our galaxy 
would produce on the order of 10,000 interactions in- 
side the Super-Kamiokande tank, yielding a rich sample 
of events for analysis. In the period of Super-K run- 
ning there has been no observed galactic supernova. In 
Ref. a limit on the rate of core-collapse supernovae 
out to 100 kiloparsecs using the SK-I and SK-II data sets 
was reported to be < 0.32/year. Although no core col- 
lapse was detected by Super-Kamiokande, it is believed 
that the universe should be bathed in the relic neutrinos 
of all of the supernovas that have exploded in the past. 
In Ref. [4(| the SK-I data sample was examined for such 
events. There is a region of energy between the endpoint 
of solar the solar neutrino spectrum and before the de- 
cay electrons of cosmic ray muons where one can hope to 
see the small number of expected events. No excess over 
background was observed, setting a limit just above the 
expected signal from many models. 

The upward-going muons used for the atmospheric os- 
cillation analysis can also be used to search for high en- 
ergy sources of neutrinos. The only muons that come 
from below the tank are generated by neutrinos, and the 
neutrinos that produce them range from approximately 
100 GeV to 100 TeV. 
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FIG. 12: Super-K 90 % CL exclusion region in WIMP param- 
eter space for a WIMP with spin-dependent coupling along 
with corresponding 90 % CL exclusion limits from UKDMC 
(dashed) and ELEGANT (dot-dashed). Figure and caption 
from Ref. El- 
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In addition to looking for neutrinos from astrophysical 
objects there is an expected neutrino signal from WIMP 
dark matter. Since WIMPs undergo gravitational attrac- 
tion they are expected to gravitate around heavy bodies 
such as the Earth and the Sun. If a WIMP and anti- 
WIMP annihilate near the center of one of these objects 
where they have accumulated, they will decay into stan- 
dard model particles some of which decay into neutrinos. 

Therefore, if there is an excess signal of upward go- 
ing neutrinos coming from either the center of the Earth 
or the Sun it can be interpreted as a signal for dark 
matter. It is shown in Ref. [4l[ that the limit on spin- 
dependent couplings inferred from the lack of an excess 
in the direction of the sun is on the order of 100 times 
more sensitive than comparable terrestrial direct dark 
matter experiments. Fig. [T^] taken from Ref. 41] shows 
a comparison with the Super-K limits with that of other 
direct-detection experiments. 

Additionally, several pure astronomical analyses have 
been performed with Super-K using samples of both 
downward-going and upward-going muons. In Ref. [42j 
the anisotropy of the primary cosmic ray flux with a mean 
energy of 10 TeV was measured by studying the relative 
sidereal variation of downward-going muons in the Super- 
K tank. 

By using upward- going muons, searches were per- 
formed for sources of high energy astrophysical neutrinos. 
In Ref. (43j the sample of upward-going muons were used 
to search the sky for point sources and signatures of a dif- 
fuse flux of neutrinos from the galactic plane. Additional 
time correlated searches with some known sources were 
also performed. Above 1 TeV of muon energy, the flux 
of upward going muons from atmospheric neutrinos be- 
comes small enough to look for a diffuse flux of extremely 
high energy neutrinos from astrophysical sources such as 
active galactic nuclei. In Ref. [44J the outer detector of 
Super-K was used to measure the direction of extremely 
high energy events. One event was found which was con- 
sistent with the background expectation and limits were 
set on this flux. 



IV. CONCLUSIONS 

This review of the published work of the Super- 
Kamiokande collaboration has shown the depth and im- 
portance of the work that the experiment has achieved. 
From astronomy and astrophysics, to tests of grand uni- 
fied models, and to the observation of neutrino mass us- 
ing neutrinos from the atmosphere and the Sun, Super-K 
has has a major impact on particle physics. 

It would be a mistake to take this snapshot of the 
published work of Super-Kamiokande as its final word. 
There is a large data set from SK-II which is currently 
being analyzed, and will soon be published. These anal- 
yses will increase the precision of the parameters Super- 
Kamiokande has already measured. In addition, the 
larger data set will allow analyses for more subtle effects 



which have not yet been studied. 

SK-I + SK-II 
6 i — ■ — ■ — ■ — i — ■ — ■ — 



5 - 




1 ' ' ' 1 ' ' ' 1 ' ' ' 1 

0.7 0.8 0.9 1 

sin 2 2e 

FIG. 13: Preliminary allowed region of atmospheric oscilla- 
tion parameters using the combined data sets of SK-I and 
SK-II. The fit to the data employed 380 bins, 70 systematic 
errors and had a Chi-squared probability of 18%. Note the 
use of the linear y-scale in the figure. 

As one example of this, Fig. Q2] shows the latest re- 
sult of the atmospheric oscillation analysis using the com- 
bined SK-I and SK-II data sets. In addition to the in- 
crease in statistics this analysis has improvements in the 
binning and systematic errors. It should be noted that 
in this figure the y-axis is plotted on a linear scale. 

This analysis employed 380 bins and 70 systematic er- 
rors. The best fit to the atmospheric oscillation parame- 
ters were: 

Am 2 = 2.5 x 10~ 3 eV 2 sin 2 29 = 1.00, (6) 

with the allowed region at 90% CL being equal to: 

1.9 x 10~ 3 eV 2 < Am 2 < 3.1 x 10~ 3 eV 2 sin 2 29 < 0.93. 

(7) 

Since June of 2006 SK-III has been operational and 
taking data. Fig [14] is a picture of SK-III during filling 
before operations commenced. Super-Kamiokande will 
continue to collect data, wait for supernovae, and in 2009 
will once again become a target for an accelerator based 
neutrino beam. 
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FIG. 14: Picture of SK-III during filling and before the beginning of operations. Photo-credit: Kamioka Observatory, 
1 CRR( Institute for Cosmic Ray Research), The University of Tokyo. 
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